Vortex Charges in High Temperature Superconductors 
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Based on an effective model Hamiltonian with competing antiferromagnetic (AF) and d-wave 
superconductivity (DSC) interactions, the vortex charge in high Tc superconductors is investigated 
by solving self-consistently the Bogoliubov-de Gennes equations. We found that the vortex charge is 
always negative when a sufficient strength of AF order is induced inside the vortex core, otherwise, 
the vortex charge is positive. By tuning the on-site Coulomb repulsion U or the doping parameter 
S, a transition between the positive and negative vortex charges may occur. The vortex charge at 
optimal doping has also been studied as a function of magnetic field. Recent NMR and Hall effect 
experiments may be understood in terms of the present results. New imaging experiments should 
be able to probe the vortex charge directly. 
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The vortex structure in high temperature supercon- 
ductors (HTS) has attracted significant interest for many 
years. Since the parent compounds are antiferromagnetic 
(AF) Mott insulators, novel physical properties of HTS 
including those in the vortex state would be expected due 
to the competition between spin magnetism and super- 
conductivity in these systems. It has been shown theoret- 
ically [0, |, |, g |, |, 0, |] that the AF order may appear 
and coexist with the underlying vortices. In a neutron 
scattering experiment by Lake et al. a remarkable 
AF-like spin density wave (SDW) was observed in the op- 
timally doped La2-a;Sra;Cu04 in the presence of a strong 
magnetic field. A moun spin rotation measurement by 
Miller et al [|lo| studied the internal magnetic field distri- 
bution in the vortex cores of underdoped YBa2Cu306+x7 
and it was revealed a feature in the high-field tail which 
fits well to a model with static alternating magnetic field. 
A very recent nuclear magnetic resonance (NMR) exper- 
iment by Mitrovic et al. ||ll|] showed that the presence 
of AF order is markedly enhanced in the vortex cores of 
near-optimally doped YBa2Cu307_5. These experiments 
have provided a strong support for the existence of AF 
order inside the vortex core in appropriately doped HTS. 

On the other hand, the vortex charge in supercon- 
ductors has also been paid considerable attention both 
theoretically ^ |l5|, |l^ and experimentally 0, 
|l^, K^. In the framework of the BCS theory, Blatter et 
al pointed out that for 5-wave superconductor the 
vortex charge is proportional to the slope of the den- 
sity of states at the Fermi level. Hayashi et al. [Q pro- 
posed that the vortex charge is always hole-like and is 
determined by the quasiparticle structure which is inde- 
pendent of the slope of the density of states. However, 
the NMR and nuclear quadrupole resonance (NQR) mea- 
surements on YBCO [|l^ seemed to obtain results for the 
vortex charge, contradictory to that predicted from the 
existing BCS theory regarding to both sign and order 
of magnitude. In view of this significant deviation, to- 
gether with the fact that the strong electron correlation 
with the d-wdive superconducting (DSC) pairing has not 



been considered in the existing theories for the vortex 
charge, we believe that the vortex charge in HTS should 
be strongly influenced by the competition effect from the 
AF and DSC orders, of which the former will play a cru- 
cial role in determining the charge nature. Also inter- 
estingly, Hall effect experiments [|l^ for HTS seemed to 
indicate that the Hall signal is electron-like in the under- 
doped up to slightly overdoped regime but hole-like in 
the overdoped regime, which could be related to the sign 
of vortex charge jl^]. Therefore, it is important to de- 
velop a sound theory for the vortex charge with the strong 
electron correlation and the (i-wave feature of HTS being 
taken into account. 

In this Letter, we shall answer two crucial questions in 
detail: what is mainly responsible for the vortex charge 
in the HTS? and how is the sign of vortex charge af- 
fected by the doping and the on-site Coulomb repulsion 
U7 Based on a widely adopted effective model Hamilto- 
nian with competing SDW and DSC orders and using a 
well-developed numerical method [|oj , we study the vor- 
tex charge in the mixed state of (i-wave HTS subjected to 
a strong magnetic field. It is found that the vortex charge 
is mainly determined by the competition of the AF order 
and the DSC order at the vortex core, and the electronic 
structure of vortex core can contain either AF order or 
normal state, corresponding to negative (electron-like) 
or positive (hole-like) charge. By tuning U or the dop- 
ing parameter, the transition between these two kinds of 
vortices occurs. 

Let us begin with an effective model Hamiltonian in a 
two-dimensional (2D)lattice, in which both the DSC and 
SDW orders are taken into account: 

+E(Aii4^su + /i.c.), (1) 

where c]^ is the electron creation operator, /i is the chem- 
ical potential, and the summation is over the nearest 



2 



neighboring sites. In the presence of magnetic field B 
perpendicular to the plane, the hopping integral can be 
expressed as tij = toex.-p[i-^ J^' A(r) • dr] for the near- 
est neighboring sites (i, j), with $o = h/2e as the su- 
perconducting flux quantum. In the presence of a strong 
magnetic field, we assume the applied magnetic field to 
be uniform and choose a Landau gauge A = {—By^ 0, 0). 
Since the internal magnetic field induced by the super- 
current around the vortex core is so small comparing 
with the external magnetic field that the above assump- 
tion is justified. The two possible SDW and DSC orders 
in cuprates are defined as Af^^ = U{cl^Ci^ — 
and Aij = VDSc{ci]Cji - Ci|Cj|)/2, where U and Vdsc 
represent respectively the interaction strengths for two 
orders. The mean- field Hamiltonian (1) can be diago- 
nalized by solving the resulting Bogoliubov-de Gennes 
equations self-consistently 
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where the single particle Hamiltonian Hij^a '- 
{Uriia- - iJi)5ij, and n^^ = Z;^K|P/(£^n) 
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(2^t)' with f{E) as the Fermi distribu- 
tion function and the electron density Ui = n^^ + n^^. 
The DSC order parameter at the ith site is Af^ = 
(Af+,^,, + Af_,^,, - Af^^^^^ - Af^^_J/4 where Ag = 

Aij exp[i^ . ^j.] a]2d Gx^y denotes the unit 

vector along (x, y) direction. The main procedure of self- 
consistent calculation is summarized as follows. For a 
given initial set of parameters Uia- and A^j, the Hamil- 
tonian is numerically diagonalized and the electron wave 
functions obtained are used to calculate the new param- 
eters for the next iteration step. The calculation is re- 
peated until the relative difference of order parameter be- 
tween two consecutive iteration steps is less than 10~^. 
The solutions corresponding to various doping concentra- 
tions can be obtained by varying the chemical potential. 

In our calculation, the length and energy are measured 
in units of the lattice constant a and the hopping integral 
to respectively. Magnetic unit cells are introduced where 
each unit cell accommodates two superconducting flux 
quanta. The related parameters are chosen as : the DSC 
coupling strength is Vdsc = 1-2, the hnear dimension of 
the unit cell of the vortex lattice is Nx x Ny = 40 x 20. 
This choice corresponds the magnetic field B ^ 37T. The 
calculation is performed in very low temperature regime. 

Our numerical results indeed show that the AF order 
is absent inside the vortex core for small U and is induced 
when U becomes larger. In Fig. 1, we plot typically the 
spatial profiles of the vortex structure for two types of 
vortices: a normal d-wave vortex core for small U (= 2.0), 
where the AF order is absent, and an AF core for larger 
U{= 2.4), where the AF order is nucleated and spreads 
out from the core center. They are obtained at the op- 
timal doping S = 0.15. Panels (a)-(c) correspond to the 



normal core while (d)-(f) for the AF core. Panels (a) and 
(d) in Fig. 1 illustrate the DSC order parameter pattern, 
which vanishes at the vortex core center. The center of 
the vortex core is situated at site (10,10). Panels (b) 
and (e) display the spatial distribution of the staggered 
magnetization of the induced AF-like SDW order defined 
as Mf = {-lyAf^^/U . No AF order is seen in the 
normal core (for U = 2.0) while the AF order exists both 
inside and outside the core (for U = 2.4) and behaves 
like a two dimensional SDW with the same wavelength 
in the x and y directions. The size of the AF core here 
is slightly enlarged than that of the normal core. The 
induced SDW order reaches its maximum value at the 
vortex core center and its spatial profile retains the same 
fourfold symmetry as that of the pure DSC case. The 
orders of DSC and SDW coexist throughout the whole 
sample. The appearance of the SDW order around the 
vortex cores strongly enhances the net electron density 
(or depletion of the hole density) at the vortex core as 
shown in panel (f). An intuitive physical understand- 
ing of positive charge for the normal vortex core can be 



D-wave Superconductivity 




Staggered IVIagnetization 




X Y 

Electron Density 




FIG. 1: Spatial variations of the DSC order parameter Af 
[(a) and (d)], staggered magnetization M/ [(b) and (e)], and 
net electron density m [(c) and (f)] in a 20 x 20 lattice. The 
left panels [(a), (b), and (c)] and the right panels [(d), (e), and 
(f)] are for [/ = 2.0 and U = 2.4, respectively. The averaged 
electron density is fixed at n — 0.85. 
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given as follows: for a particle-hole asymmetric system 
like doped cuprates, the chemical potential for electrons 
in DSC state would be slightly lower than that of the nor- 
mal state, when a normal core is imbedded into the DSC 
background and in order to reach equilibrium, electrons 
have to flow from the inside to the outside of the core 
which leads to the electron depletion inside the vortex 
core, as shown in panel (c); while in the case of AF core, 
the hole number is suppressed and as a result, the vor- 
tex carries negative charge. The enhancement of electron 
number inside the AF vortex core has also been numeri- 
cally obtained by other calculations ||, 0, D • 

To examine the vortex charge Qy as functions of both 
5 and /7, the upper right inset in Fig. 2 plots the phase 
diagram of 5 versus U for positively (hole-rich) and nega- 
tively (electron-rich) charged vortices. It is obvious that 
the AF vortex core can easily show up in the underdoped 
regime or with stronger AF interaction while normal core 
tends to exist in the overdoped regime or with weaker AF 
interaction. The electron density inside the core is higher 
than the average density in the underdoped region while 
the electron density becomes lower than the average in 
the overdoped region. There exists a clear boundary be- 
tween these two phases. The AF order is generated in the 
region where the DSC order parameter is suppressed. To 
estimate the core charge of a single vortex, we first deter- 
mine the vortex size by examining the spatial profile of 
DSC order parameter. Next we make a summation of the 
net electron density inside the vortex core. As shown in 
Fig. 2, the ^-dependence of Qv/e (the electron number) 
for [/ = 2.4 exhibits a first-order like transition at ^ = (5c 
(~ 0.18). The magnitude of the discontinuity reduces to 
one third when U = 2.2. The critical value of the doping 
level 5c is [/-value dependent or the sample-dependent. 
The larger U case corresponds to larger 5^ Recent NMR 
experiments indicated that the AF order exists in 
the vortex core at the optimal doping level in cuprates, 
which may imply that the critical doping level 5c could 
be extended to slightly overdoped region. Therefore, It 
is clear that the related phenomena in the slightly over- 
doped sample may be qualitatively the same as those in 
the underdoped sample, e.g., the slightly overdoped sam- 
ple has the electron-rich vortex core as well. This result 
agrees well with the experiment for slightly overdoped 
YBa2Cu307 [|l9|, in sharp contrast to the hole-rich vor- 
tex core predicted by the BCS theory. Also interestingly, 
even though the origin of Hall sign anomaly is still de- 
batable , the vortex charge could make an additional 
contribution to the sign change in the mixed state Hall 
conductivity Our calculations which is schemati- 

cally shown in the lower left inset of Fig. 2 would favor 
that the Hall signal is electron-like from the underdoped 
to slightly overdoped regime but hole-like in the appre- 
ciable overdoped regime. This result is consistent with 
the phase diagram obtained by the Hall effect measure- 
ments 0. 

In addition, the charge magnitude estimated from the 
BCS theory pi] is two orders smaller than that of exper- 




FIG. 2: Doping dependence of the number of vortex charge 
Qv/e for U — 2A where the electron charge e < 0. The left 
inset shows the doping dependence of the sign of the vortex 
charge ( positive for hole-rich and negative for electron-rich). 
5°^^ and 5c denote respectively the optimal doping and crit- 
ical doping. The right inset represents the phase diagram of 
doping level versus interaction strength U for positive and 
negative charged vortex. 
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FIG. 3: The interaction strength dependence of the staggered 
magnetization Mc (solid line) at the vortex center and the 
number of vortex charge Qv/e (dotted line). The averaged 
electron density is fixed at n = 0.85. The inset shows the 
extrapolation of Qv/g and with magnetic field strength 
in unit of Tesla. 



imental observation for HTS. The magnitude of vortex 
charge estimated from our calculation is about 0.06e at 22 
Tesla, which seems much larger than the experimental es- 
timation 0.005e - 0.02e at 9.4 Tesla for YBa2Cu307 [0. 
The reason appears to be mainly due to a much higher 
magnetic field used in our calculation, which will lead 
to a larger AF order. The inset of Fig. 3 represents 
the approximate extrapolation of the vortex charge mag- 
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nitude versus the magnetic field. The estimated vortex 
charge at 9.4 Tesla is indeed in the same order of mag- 
nitude as reported in the experiment. From Fig. 3, one 
can clearly see an abrupt jump for the number of vor- 
tex charge Q^/e and staggered magnetization at the vor- 
tex core center as U varies around 2.11, and this 
positively-negatively charged vortex transition appears 
also to be first-order like. It is now quite clear that the 
vortex charge is strongly influenced by two competing 
effects — the suppression of the DSC order at the core 
center which leads to the depletion of the electrons and 
the induction of the AF order which favors the accumu- 
lation of electrons. Whether the negative vortex charge 
appears depends solely on whether there is a sufficient 
AF order inside the vortex core, as clearly seen in Fig. 3. 
Although our calculation is based upon the phenomeno- 
logical Hamiltonian, our results are robust despite of dif- 
ferent band parameters and should give a qualitative de- 
scription on the vortex physics in HTS. 

We now turn to discuss the experimental results for 



strongly underdoped YBa2Cu408 |19 where a positively 
charged vortex is reported. This seems to be inconsistent 
with our prediction. We believe that the vortex charge in 
the above experiment was deduced from an oversimplified 
assumption that the electron density is uniform either in 
the absence of magnetic field or far away from the vortex 
core in this strongly underdoped HTS. In fact, experi- 
ments showed clearly the remarkable inhomogeneities in 
the underdoped sample |2^, Many theoretical stud- 
ies including the present one also show the presence of 
stripe-like charge density structures in the strongly un- 
derdoped sample ||, |^, U, |5|. Upon the application 
of a magnetic field, the spatial charge distribution could 
become more inhomogeneous even away from the vortex 



core. Therefore, their estimation of the vortex charge 
for the underdoped YBa2Cu408 might be invalid. For 
the slightly overdoped HTS, the sample is less inhomoge- 
neous, and their estimation may be qualitatively correct. 

With respect to the complexity of the underdoped case, 
we suggest to use the spatially resolved high magnetic 
field NMR to probe the vortex charge. In this way, a 
clear resolution of the vortex core region can be reached 
since the fraction of the spectrum inside the core grows 
with the increase of the magnetic field. It seems better 
to probe the vortex charge in slightly underdoped sam- 
ples to test our results because the strongly underdoped 
samples have the complications mentioned above. We 
would also like to pinpoint that the high resolution STM 
may be a good candidate to probe the vortex charge by 
integrating local density of states up to the chemical po- 
tential. The spatial electron density distribution can also 
be directly imaged by the electrostatic force microscope, 
which detects the force gradient acting on the tip, and 
the scanning surface potential microscopy, which mea- 
sures the first harmonic of the force. If the vortex indeed 
possesses the charge as we find here, these direct imaging 
techniques can be utilized as powerful tools to study the 
vortex dynamics in HTS. 
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